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_ , Abstract 

m 

' The spectrum of B and D mesons (including the low lying orbitally and radially ex- 

. cited states) is calculated using the quark-antiquark Hamiltonian derived from QCD in 

K*- i the einbein field formalism. Spin-spin and spin-orbit terms due to the confinement and 

' OGE interactions are taken into account as perturbations. Results for the masses and 

■ splittings are confronted to the experimental and recent lattice data and are demon- 
OO , strated to be in a reasonable agreement with both. We find that the orbital excitations 

with / = 2 and I = 3 for D meson lie approximately in the same region as its first 

■ radial excitation that might solve the mystery of the extremely narrow D(2637) state 
, recently claimed by DELPHI Collaboration. 



&: 1 Introduction 



In spite of a rather long history of discussion and many theoretical attacks, an interest to the 
^ . properties of heavy-light mesons is still very high. Indeed, heavy-light systems incorporate 
both properties of the light quarks which are extremely important for the physics of Chiral 
Symmetry Breaking (CSB), one of the most challenging phenomenon of the nonperturbative 
QCD, and heavy degrees of freedom which allow application of such profound methods of 
investigation as Heavy Quark Effective Theory (HQET) or Operator Product Expansion 
(OPE). The two main sources of information, experiment and lattice simulations, deliver 
new data on the heavy-light mesonic spectra, including orbital and radial excitations of 
the qq pair, which strongly need theoretical identification and description. In the present 
paper we calculate the masses of several low lying D and B mesonic states using the method 
developed in g based on the Hamiltonian approach to the bound states problem in QCD. 
Starting with the quarks with current masses we intensively use the so-called einbein fields 
formalism P, ^ which not only considerably simplifies calculations, but allows to justify the 
treatment of the spin-spin and spin-orbit interactions as perturbations, naturally explaining 
the appearance of the "constituent" quark mass. Results for the spectra and splittings are 
compared with the experimental and lattice data and a good agreement is found. 

*e-niail: yulia@heron.itep.ru 
^e-mail: nefediev@heron.itep.ru 
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We demonstrate that in our model the proper account of the QCD string dynamics 
gives extra negative contribution into the energy of the orbital excitations and makes the 
orbitally excited state (/ = 2, n = and / = 3, n = 0) share the region of masses with the 
first radial excitation (/ = 0, n = 1). The latter observation may resolve the problem of 
the resonance D(2637) recently claimed by DELPHI Collaboration ^ and which is under 
intensive discussion at present (see e.g. 0). Indeed, in spite of the fact that in our model both 
above mentioned states lie somewhat higher than the experimentally observed resonance, the 
latter can be identified with the orbital excitation of charm (with / = 2 or / = 3) rather then 
with the first radial excitation (n = 1) solving in such a way the problem of an extremely 
small width of the resonance inconsistent with any estimates of the 2S'-state decays [0. 



2 Einbein fields formalism and the mesonic Hamilto- 
nian 

The einbein fields formalism which has a long history in literature was introduced as a 
method allowing to treat the kinematics of relativistic particles, but then it was generalized 
for the case of spinning particles and strings p, ^ . One of its applications is the possibility 
to develop the canonical quantization procedure a la Dirac |^. In its simplest form this 
formalism can be applied to the point-like scalar relativistic particle described by the action 

S = f Ldt L = —m\/~^ = (1) 

so that the modified form of (0) looks like 

where the einbein field /i is introduced. Dynamics defined by the equations of motion for 
Lagrangians (|l]) and (0) is the same if the einbein field is treated as an independent degree 
of freedom and the corresponding constraints are introduced and properly treated. For a 
short review of the einbein fields formalism see e.g. and references therein. In the mean 
time there exists another approach to the einbein fields 0| which allows to neglect their 
dependence on the proper time t and thus to treat them as variational parameters to be 
got rid of in the spectrum rather than in the Hamiltonian. Such an approach was used to 



a success in a set of papers (see e.g. |TT|) and proved to be rather accurate reproducing 
the relativistic spectra with the error within several per cent. In what follows we accept the 
above accuracy and treat einbeins variationally. 

As it was shown in p, |^ within the Vacuum Correlators Method (VCM) [0, writing 
the gauge invariant Euclidean Green's function of the qq meson 

Ggg = {^+^{x, (x, y\A))A ^gg{xuX2\A) = ^5(a;i)$(xi, a;2)^,(x2), (3) 

with $(xi,X2) being the standard path-ordered parallel transporter operator 

^{xi,X2) = Pexp (^ig J dZf^A^,^ 

and employing Feynman-Schwinger representation for the single particle Green's functions 
and the area law for the isolated Wilson loop bounded by the quark and antiquark trajec- 
tories, one can arrive at the following result P] 



Ggg = J Dfi,{ti)Dfi2{t2)DxiDx2e-'''-'''Tr 



r(^) (mi - Z))r» (m2 -b)x (4) 
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P ( (1) ^ ^ / _dt2_ (2) ^ ^ 

where fj,i and /i2 are the einbein fields, Ki being the kinetic energies of the quarks 

A'.= r*f?^ + ^ + ^V >=1,2, (5) 




S" stands for the minimal area swept by the quark- ant iquark trajectories, cr^,^ = j-i'^^'yu — 
Ivlpi) and 5/5Sfj_y denotes the derivative with respect to the element of the area S. The 
form of the kinetic energy (^) immediately follows from the equation (|]) if the laboratory 
gauge t = a^o is fixed in the proper time reparametrization group. Equation allows to 
find the Hamiltonian of the quark-antiquark pair connected with the relativistic string. If 
one neglects spin-dependent terms in (^, then the corresponding Hamiltonian coincides by 
the one for the Nambu-Goto string with scalar quarks at the ends which follows from the 
Lagrangian 

L{t) = —mi^J±l — rn2\/i!l — a dl3\J {ww'Y — w'^w''^ xio = X2o = t, (6) 



where the minimal string is parametrized by the profile function w{t,j3), the dot and the 
prime denote its derivatives with respect to the proper time t and the string coordinate /3 
respectively. Thus for the centre-of-mass frame one arrives at 




2^. 2) Jo ^V2^ 2; 2r2[/xi(l-C)' + /i2C' + /o W/3-C)T 

A^i + Ai2 + Jo dPu 

Note that an extra einbein u following equation (|^ has the meaning of the string energy 
density. Extrema in all three einbeins {fii, fi2 and u) should be taken in (0). For the case of 
zero angular momentum L = the last term in (|^) vanishes and the extremal value for u is 
easily found to be 

z/Q = err, (8) 

so that the resulting interaction gives the linearly rising potential ar. Expansion of the last 
term in Hamiltonian in powers of \fo j with the substitution of the extremal value z/q 
yields 

H = Ho + Vstring (9) 

-o^K^^D-^-i^-Co, (10) 

Vstring g^2^2 ^ 

where we have supplied Hq with the constant term Co and the colour Coulomb interaction 
which comes from the perturbative gluon exchange omitted in (|^) for simplicity. Equation 
( |TID gives the first spin-independent correction called the string correction |jT3|, |l| which is 
known to reduce the inverse Regge trajectory slope dMl/dL = 4cr specific for the heavy-light 
limit of the Hamiltonian Hq, and, when taken at the account at full scale, to lead to the 



correct string slope na [||, [TJ . 



3 



In what follows we shall consider Ho as the zeroth approximation for calculation of the 
mesonic spectrum. Once the angular momentum L, the total spin of the quark-antiquark 
pair S and the total momentum J = L + S are separately conserved with the Hamiltonian 
Hq then one can specify the mesonic states as terms where n is the radial excitation 

number. 

Let us turn back to the equation and to identify the spin-dependent correction to the 
zeroth order Hamiltonian Hq. Following [Q] one finds 

2 

+ -j^^ (3(5in)(52n) - (8182)) + {SL) (2 - ln(/ir) - 7^), = 0.57 (12) 

where k = and we have added the term of order which comes from one-loop cal- 




culations and is extensively discussed in literature M, llH]. Note that up to the last term 



potential (|T^) coincides in form with the Eichten-Feinberg-Gromes results [T^, but con 



tains /i's in the denominators instead of masses. In the meantime as clearly seen from (|1^ 



/^i ~ Y ii^) + mf > rrii or even yUj ^ rrij for light quark^ (see also Table 1). Einbeins /I's 
can be viewed as dynamical or "constituent" quark masses that makes it possible to consider 
(|12|) as a correction and to justify the expansion made in (|1T]) (see e.g. for more detailed 
discussion of this issue). 



3 Numerical solution 



3.1 Spectrum of the Hamiltonian Hq 

As mentioned in the introduction einbein fields can considerably simplify the relativistic 
dynamics. Indeed with the einbeins /I's introduced the kinetic part of the Hamiltonian Hq 
from ([T0|) has a nonrelativistic form, so that the corresponding Schrodinger equation can be 
written in the reduced dimensionless form (see e.g. [jl8|, |17| ) 



1^1 ^\ ^ 2/i \ 

+ FI - 1^ I Xa = a{\)xx A = K I I 



2/3 



(13) 



On finding the solutions of equation ([13|) for x\ and a(A) as functions of the reduced 
Coulomb interaction strength A one has for the extremal values of the einbeins: 



/ii(A) = v'mf + A2(A) /i2(A) = V'mi + A2(A) /x(A) 



1 



3/2 



(14) 



where 



A^(A) 



3k 



a + 2X 



da 



dX 



Then the definition of the reduced field fi gives the equation for A 

/ii(A)/i2(A) 



/i(A) 



/^i(A)+/^2(A)' 



(15) 



^Let us remind here that these are current quark masses to be denoted m.i and to enter equations (p])-(pl] 
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n 


/ 


meson 


nil 


m2 


a 


as 


A 


/ii 


1^2 


/i 




|V^(o)| 








D 


1.4 


0.009 


0.17 


0.4 


0.817 


1.497 


0.529 


0.391 


2.198 


0.161 






Ds 


1.4 


0.17 


0.17 


0.4 


0.847 


1.501 


0.569 


0.412 


2.224 


0.167 






B 


4.8 


0.005 


0.17 


0.39 


0.999 


4.840 


0.619 


0.549 


5.527 


0.209 






Bs 


4.8 


0.17 


0.17 


0.39 


1.035 


4.842 


0.658 


0.579 


5.550 


0.219 





1 


D 


1.4 


0.009 


0.17 


0.4 


0.869 


1.522 


0.597 


0.428 


2.640 









Ds 


1.4 


0.17 


0.17 


0.4 


0.891 


1.525 


0.629 


0.445 


2.663 









B 


4.8 


0.005 


0.17 


0.39 


1.052 


4.847 


0.675 


0.593 


5.949 









Bs 


4.8 


0.17 


0.17 


0.39 


1.080 


4.849 


0.707 


0.617 


5.970 
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D 


1.4 


0.009 


0.17 


0.4 


0.924 


1.554 


0.674 


0.470 


2.961 









Ds 


1.4 


0.17 


0.17 


0.4 


0.942 


1.557 


0.702 


0.484 


2.982 









B 


4.8 


0.005 


0.17 


0.39 


1.128 


4.860 


0.762 


0.659 


6.245 









Bs 


4.8 


0.17 


0.17 


0.39 


1.151 


4.861 


0.789 


0.679 


6.263 





1 





D 


1.4 


0.009 


0.17 


0.4 


0.929 


1.557 


0.682 


0.474 


2.848 


0.162 






Ds 


1.4 


0.17 


0.17 


0.4 


0.947 


1.561 


0.710 


0.488 


2.869 


0.165 






B 


4.8 


0.005 


0.17 


0.39 


1.142 


4.863 


0.779 


0.671 


6.131 


0.207 






Bs 


4.8 


0.17 


0.17 


0.39 


1.165 


4.864 


0.806 


0.692 


6.149 


0.212 



Table 1: Solutions of the equations (P!^)-([T3|) for standard values of the string tension a, the 
strong coupling constant and the masses of the quarks. All parameters are given in GeV 
to the appropriate powers. 



which is the subject to numerical investigation. Unfortunately there is no analytic solution 



for a(A), so one has to generate a selfconsistent solution of both equations, (|l^) and dlSf) . 



Let us give here two more formulae which will be needed for calculating the spin-spin 



and spin-orbit splittings. For radially excited states one can find |I8|, |17 



1^(0)1^ = ^(l + A(x-)), (16) 

where _ 

(r^) = (2/ia)^/=^(a;^) = (2/1^)^/^ / x^^^ |^^^^)|2 ^ ^ > _ 3/. (17) 







In Table 1 we give the results of numerical calculations with the Hamiltonian Hq for D, 
Ds, B and Bs mesons. 



3.2 Spin-spin and spin-orbit splittings. Comparison with the ex- 
perimental and lattice data 

In this subsection we calculate the contribution of the Ktring and Vgd terms. First of all we 
collect the averaged values of the spin-orbit and spin-tensor interactions between various IP 
and ID eigenstates of the Hamiltonian Hq 

(iPi|5iL|iPi) = (iPi|52L|iPi) = CPi|(5in)(52n)|iPi) = -1 



3 

(3Po|5iL|3Po) = -1 (3Po|S^2L|3Po) = -1 (3p„|(Siri)(S2n)|3Po) = -| 
(3Pi|5iL|3Pi) = -1 (3Pi|S^2^|3p,) = _l {^Pi\{Sin){S2n)fPi) = l 

f P2|5iL|3P2) = 1 (3P2|S^2^|3P2) = ^ fPs | (5^1^) (5^2^) 1 3P2) = ^ 
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;i8) 



{^D2\SiL\^D2) 


= {^D2\S2L 


'D2) 


= (iD2|(Sirr)(52n) 


'D2) 


{^D,\S,L\^D,) 


= -| ('D,\S2L 




= -| (3Di|(^in)(^2fz) 




(^D2\SiL\^D2) 


= — -i- Do \S'}L 


^D2) 

^1 


= -i (3L)2l(M(M 




{^D,\SiL\^Ds) 


= 1 (3^315^2^ 




= 1 fZ)3|(Sir?)(52n) 




as well as the transition matrix elements 








CPi|S^iL|3Pi) = 


1 

V2 


CPi|s^2^pPi) = -4= 

V ^ 






{^D2\SiL\^D2) = 


/3 
' V2 


(1^2152^13^2) = -Wf 





1 

'5 



(19) 



1 

2S 



(20) 



which lead to mixing of I^Pi), \^Pi), and 1^-02); P-D2) pairs so that the physical state is 
subject to the matrix equations of the following type: 



Ei-E V12 
Vy2 E2 — E 



(21) 



The results of our numerical calculations as well as the experimental and lattice data are 
given in Tables 2,3. The only fitting parameter we use here is the constant Co which takes 
the following values: 

Co{D) = 212MeV Cq{D,) = 124MeV Co{B) = 203MeV Co{Bs) = 124MeV. (22) 



As clearly seen from ( p^ the constant Cq practically does not depend on the heavy quark 
and is completely defined by the properties of the light one. 

4 Discussion and conclusions 

Several comments concerning our choice of parameters are in order here. From Table 1 it 
is seen that we use practically the same values of the strong coupling constant for all 
states. To justify this action let us note that in the case of a light particle moving in the 
field of a heavy one the OGE interaction between them depends not on the total mass of 
the system, but rather on its size. Then the slight variation of the strong coupling is due to 
small differences in the mesons sizes. The exact value of a, is chosen to be near its frozen 



value 1 19] that emphasizes the nonperturbative nature of the processes responsible for the 
formation of the heavy-light mesons spectrum. The fits are very weakly sensible to variations 
of the heavy quarks mass which lead mainly to rescaling of the overall shift constant Cq, 
so that the latter remains practically the only fitting parameter as stated in the preceding 
section. 

In conclusion we would like to comment on the possible identification of the resonance 
D(2637) recently claimed by DELPHI Collaboration 0. It was reported to be extremely 
narrow, about 15MeV, that leads to a confusion as such a small width is inconsistent with 
the identification of this state as the first radial excitation D*\J^ = It was found in a 
number of papers that in spite of the fact that the mass of this state perfectly coincides 
with the predictions of quark models, all estimates of its width fail to give such a small value. 
In the meantime it was found that widths of orbitally excited D mesons with quantum 
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Table 2: Masses of the D, Dg, B and Bg mesons in MeV . Lattice results are extracted from 
Figures 26,27 and Tables XXVIII,XXIX of l2|. Symbols l^Pi/^Pi and l^D2/^D2 are used 
to indicate that the physical states are mixtures of the l^Pi and l^Pi or 1^D2 and 1^D2 
states correspondingly. Underlined figures give masses of the most probable candidates for 
the experimentally observed resonances. 



Splitting 


Ds-D 


D*g - D* 


D* - D 


D: - Dg 


Bs-B 


B* - B* 


B* -B 


b: - Bs 


Experim. 


99 


102 


141 


144 


90 


91 


46 


47 


Theory 


114 


115 


146 


147 


100 


102 


63 


65 


Lattice 


100 


92 


110 


103 


90 


90 


30 


29 



Table 3: Splittings for the D, Ds, B and Bs mesons in MeV. Lattice results are taken from 
Tables XXVIII,XXIX of [EDI. 
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numbers 2^, 3^ could be consistent with the reported value J^, but then the following two 
points were put forward as main objections: i) a neighboring slightly more massive state 
should be observed as well, ii) quark models predict orbitally excited mesons to be at least 
bOMeV heavier than needed. Here we would like to comment on the second argument. As 
one can find from the fit given in Table 2, in our model the first radial and the second 
orbital excitations share the same region of masses. It is not surprise as the negative mass 
shift of about 50MeV for the orbitally excited state, missing in standard quark models, is 
readily delivered by the string correction given by equation (pAl). This contribution does not 
affect radially excited states but is significant for orbital excitations. It was demonstrated 
in that the account of the proper dynamics of the QCD string in mesons brings the 



slope of Regge trajectories to their correct values. Appearance of the term (|TT]) in the meson 
Hamiltonian and its important contribution into the masses of orbitally excited states is yet 
another reflection of the general situation that the proper QCD string dynamics is extremely 
important in description of hadronic properties and thus it should be taken into account. 
So we conclude that the mysterious D(2637) state indeed can be identified with the orbital 
excitation 2~ or 3~ rather than with the radial one 0~, that resolves the second objective 
above. This statement may hold true for the corresponding states in the B meson spectrum, 
where DELPHI Collaboration also claims a similar state ||2T[| . 

As far as the first one is concerned, our model also predicts another orbitally excited state 
with the mass 2728Mey, i.e. 65MeV higher than D2 given in the table, and its experimental 
grounds are really not clear. 

Our predictions for D(2654), D(2663) and Z)(2664) lie somewhat higher than the exper- 
imentally observed value. We find it to be a reflection of a general lack of the "/z-technique" 
used in this paper, which gives larger errors for higher excited states. The difference of about 
15 — 20MeV between the theoretical predictions and the experimental datum 2637± 6Mel^ 
can be explained in this way. Development of a systematical approach to the einbein fields 
as variational parameters could shed light on the sources of systematical errors and possibly 
to improve the results presented in this paper. 

The authors would like to thank Yu.A.Simonov for constant interest to the work and 
useful discussions, R.N.Faustov for reading the manuscript and for valuable comments, and 
B.L.G.Bakker for providing the numerical solutions for the reduced Schrodinger equation. 
Financial support of RFFI grants 00-02-17836 and 00-15-96786 and INTAS-RFFI grant IR- 
97-232 is gratefully acknowledged. 



References 

[1] A.Yu.Dubin, A.B.Kaidalov, Yu.A.Simonov, Phys.Lett. B323 (1994) 41; Phys.Lett. 
B343 (1995) 310; Nuovo Cim. 107A (1994) 2499; Phys.Atom.Nucl. 58 

[2] Yu.A.Simonov, Nucl.Phys. B307 (1988) 512; 

Yu.A.Simonov and J.A.Tjon, Ann.Phys. (NY) 228 (1993) 1 

[3] L.Brink, P.Di Vecchia, P.Howe, Nucl.Phys. B118 (1977) 76 

[4] Yu.A.Simonov, Phys.Lett. B226 (1989) 151 

[5] Yu.S.Kalashnikova and A.V.Nefediev, Phys.Atom.Nucl. 60 (1997) 1389 
[6] DELPHI Collaboration, P.Abreu et.al, Phys.Lett. B426 (1998) 231 



8 



[7] D.Melikhov and O.Pene, Phys.Lett. B446 (1999) 336; 
P.R.Page, Phys.Rev. D60 (1999) 057501 

A.M.Polyakov, "Gauge Fields and Strings", Harwood Academic Publishers, 1987; 

E.S.Fradkin and D.M.Gitman, Phys.Rev. D44 (1991) 3230; 

G. V.Grygoryan and R.P.Grigoryan, Phys.Atom.Nucl. 53 (1991) 1737; 
Yu.S.Kalashnikova and A.V.Nefediev, Phys.Atom.Nucl. 62 (1999) 377 

P.A.M.Dirac, "Letures on Quantum Mechanics", Belter Graduate School of Science, 
Yeshiva University, New York (1964) 

V.L.Morgunov, A.V.Nefediev and Yu.A.Simonov, Phys.Lett. B459 (1999) 653 

H. G.Dosch, Phys.Lett. B190 (1987) 177; 
H.G.Dosch and Yu.A.Simonov, Phys.Lett. B205 (1988) 339; 
Yu.A.Simonov, Nucl.Phys. B307 (1989) 67 

J.Merlin and J.Paton, J.Phys. Gil (1985) 439 

J.Pantaleone, S.-H.Henry Tye, Y.J.Ng, Phys.Rev. D33 (1986) 777 

Yu.A.Simonov, Nucl.Phys. B592 (2000) 350 

E.Eichten and F.L.Feinberg, Phys.Rev. D23 (1981) 2724; 
D.Gromes, Z.Phys. C26 (1984) 401 

Yu.A.Simonov, Lectures given at the XVII International School of Physics "QCD: Per- 
turbative or Nonperturbative" , Lisbon, 29 September — 4 October 1999 

Yu.A.Simonov, Z.Phys. C53 (1992) 419 

Yu.A.Simonov, Yad.Fiz. 58 (1995) 113; JETP Lett. 57 (1993) 513 
J.Hein et.al, Phys.Rev. D62 (2000) 074503 

M.Fcindt, O.Podorin (DELPHI Collaboration), contr.paO 1-021 ICHEP'96 Warshaw, 
DELPHI 96-93 CONE 22; 

C.Weiser, Proceedings of the 28th International Conference on Hight Energy Physics, 
editors Z.Ajduk, A.K.Worblewski, (World Scientific, Singapore 1997), page 531 



9 



